Draft version February 5, 2008 

Preprint typeset using IATj^X style cmulateapj v. 10/09/06 



& K. KUCHIBHOTLA 



o 
o 

(N 

a 

i> 

(N 

S 1 
6 



> 
in 



O 

o 

> 

x 



TMC-1C: AN ACCRETING STARLESS CORE 
S. Schnee 1-2 , P. Caselli 1,3 , A. Goodman 1 , H. G. Arce 4 , J. Ballesteros-Paredes 5 

Draft version February 5, 2008 

ABSTRACT 

We have mapped the starless core TMC-1C in a variety of molecular lines with the IRAM 30m 
telescope. High density tracers show clear signs of self-absorption and sub-sonic infall asymmetries 
are present in N2H + (l-0) and DCO + (2-l) lines. The inward velocity profile in N2H + (l-0) is extended 
over a region of about 7,000 AU in radius around the dust continuum peak, which is the most extended 
"infalling" region observed in a starless core with this tracer. The kinetic temperature (~ 12 K) 
measured from C 17 and C 18 suggests that their emission comes from a shell outside the colder 
interior traced by the mm continuum dust. The C ls O(2-l) excitation temperature drops from 12 K 
to ~10 K away from the center. This is consistent with a volume density drop of the gas traced by the 
C 18 lines, from ~4xl0 4 cm -3 towards the dust peak to ~6xl0 3 cm -3 at a projected distance from 
the dust peak of 80" (or 11,000 AU). The column density implied by the gas and dust show similar 
N2H + and CO depletion factors (fr> < 6). This can be explained with a simple scenario in which: 
(i) the TMC-1C core is embedded in a relatively dense environment (n(Ha) — 10 4 cm -3 ), where CO 
is mostly in the gas phase and the N2H + abundance had time to reach equilibrium values; (ii) the 
surrounding material (rich in CO and N2H + ) is accreting onto the dense core nucleus; (hi) TMC-1C is 
older than 3xl0 5 yr, to account for the observed abundance of N2H + across the core (~10~ 10 w.r.t. 
H2); and (iv) the core nucleus is either much younger (~ 10 4 yr) or "undepleted" material from the 
surrounding envelope has fallen towards it in the past 10,000 yr. 
Subject headings: stars: formation — dust, extinction — submillimeter, molecules 



1. INTRODUCTION 

Dense starless cores in nearby low-mass star-forming 
regions such as Taurus represent the simplest areas in 
which to study the initial conditions of star formation. 
The dominant component of starless cores, H2, is largely 
invisible in the quiescent interstellar medium, so as- 
tronomers typically rely on spectral line maps of trace 
molecules and continuum observations of the thermal 
emission from dust to derive their kinematics and phys- 
ical state. However, it is now well established that 
different species and transitions trace different regions 
of dense cores, so that a comprehensive multi-line ob- 
servations, together with detailed millimeter and sub- 
millimeter continuum mapping are required to under- 
stand the structure and the evolutionary status of an 
object which will eventually form a protostar and a pro- 
toplanetary system. 

Previous studies of starless cores in Taurus, as well 
as other nearby star-forming regions, have shown that 
the relative abundance of many molecules varies signifi- 
cantly between the warmer, less dense envelopes and the 
colder, denser interiors (see Ceccarelli et al. 2006 and 
Di Francesco et al. 2006 for detailed reviews on this 
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topic) . Fo r instance. ICaselli et al.l (|1999f) . iBergin et al.1 
(2002J) and lTafalla et al.l (|2004) have shown that carbon- 
bearing species such as C 17 0, C ls O, C 34 S, and CS are 
largely absent from the cores L1544, L1498 and L1517B 
at densities larger than a few 10 4 cm" 3 , while nitrogen- 
bearing species such as N2H + and ammonia are prefer- 
entially seen at high densities. The chemical variations 
within a starless core are likely the result of molecu- 
lar freeze-out onto the surfaces of dust grains at high 
densities and low temperatures, followed by gas phase 
chemical processes, which are profoundly affected by the 
abunda nce drop of important species, in particular CO 
(see e.g. iDalgarno fc Lepr] 119841: IBergin fc Langerlll997t 
iTavlor et al.lll998l : lAikawa et al.ll2005h . 

In the past few years it has also been found that not all 
starless cores show a similar pattern of molecular abun- 
dances and physical structure. Indeed, there is a sub- 
sample of starless cores (often called pre-stellar cores), 
which are particularly centrally concentrated, that shows 
kinematic and chemical features typical of evolved ob- 
jects on the verge of star formation. These features in- 
clude large values of CO depletion and deuterium frac- 
tionation, and evidence of "central" infall, i.e. pres- 
ence of infall asymmetry in high density tracers in a re- 
stricted region surrounding the mm continuum dust peak 
(I Williams et~ai . 1999; Caselli et al. 2002a; Re dman et all 
I2002t ICrapsi et all 120051 : 1 Williams et~aT1l2006h . It is in- 
teresting that not all physically evolved cores show chem - 
ical ly evolved compositions, as s hown bv lLee et al.l (|2003f ) 
and rTafalla fc Santiago! (|2004bl ). It is thus important to 
study in detail a larger number of cores to understand 
what is causing the chemical differentiation in objects 
with apparently similar physical ages. This is why we 
decided to focus our attention on TMC-1C, a dense core 
in Taurus, with physical properties quite similar to the 
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prototypical pre-stellar core L1544 (also in Taurus), to 
study possible differences and try to understand their 
nature. 

TMC-1C is a starless core in the Tau rus molecular 
cloud , with a distance estimated at 140 pc (jKenvon et al.l 
1994). In a previous study, we have shown that TMC- 
1C has a mass of 6 M Q within a radius of 0.06 pc from 
the column density peak, which is a factor of two larger 
than the virial mass derived from the N2H + (l-0) line 
width, and we have sh own that there is evidence for sub- 
sonic inward motions (|Schnee &: Goo dman 20051) as well 
as a velocity gradient consistent with solid body rota - 
tion at a rate of 0.3 km s _1 pc -1 dGoodman et al.lll993f) . 
TMC-1C is a coherent core with a roughly constant ve- 
locity dispersion, slightly h igher than the sound speed , 
over a radius of 0.1 pc (|Barranco fc Goodman! 119981 : 
iGoodman et~aTl fl99l . Using SCUBA and MAMBO 
bolometer maps of TMC-1C at 450, 850 and 1200 ^m, 
we have mapped the dust temperature and column den- 
sity and shown that the dust t emperature at the center 
of the core is very low (~6 K) (|Schnee et al.ll200l . 

In order to disentangle the physical and chemical in- 
formation that can be gleaned from a combination of 
gas and dust observations of a dense core, we have 
now mapped TMC -1C at three continuum wavelengths 
(Schn ee et al.ll2007t ) and seven molecular lines. In Sec.[H 
continuum and line observations are described. Spectra 
and maps are presented in Sec. [3J The analysis of the 
data, along with the discussion, has been divided in three 
parts: kinematics, including line width variations across 
the cloud, velocity gradients and inward velocities, is in 
Sec.[4j gas and dust column density and temperature are 
in Sect. molecular depletion and chemical processes 
are discussed in Sect. [5] A summary can be found in 
Sect. [3 

2. OBSERVATIONS 

2.1. Continuum 

To map the density and temperature structure of 
TMC-1C, we have observed thermal dust emission at 
450 and 850 fim with SCUBA and at 1200 fim with 
MAMBO-2. 

2.1.1. SCUBA 

We observ ed a 10'xlO' region around TMC-1C us- 
ing SCUBA ([Holland et al.|[i999t ) on the JCMT at 450 
and 850 ^m. We used the standard scan-mapping 
mode, recording 450 and 850 /xm data simul taneously 
(jPierce-Price et all 120001: iBianchi et all I2000D . Chop 
throws of 30", 44" and 68" were used in both the right 
ascension and declination directions. The resolution at 
450 and 850 /jm is 7.5" and 14" respectively. The ab- 
solute flux calibration is ~12% at 450 /xm and ~4% at 
850 /im. The noise in the 450 and 850 /im maps are 13 
and 9 mJy/beam, re spectively. The da ta reduction is 
described in detail in lSchnee et al.l (|2007f) . 

2.1.2. MAMBO-2 

Kauffmann et al. (in prep.) used the MAMBO-2 array 
(jKrevsa et al.l ll999) on the IRAM 30-meter telescope on 
Pico Veleta (Spain) to map TMC-1C at 1200 itm. The 
MAMBO beam size is 10"7. The source was mapped on- 
the-fly, chopping in azimuth by 60" to 70" at a rate of 2 



Hz. The absolute flux calibration is uncertain to ~10%, 
and the noise in the 1200 /im map is 3 mJy/beam. The 
data reduction is described in detail in Kauffmann et al. 
(in prep.). 

2.2. Spectral Line 

We have used the IRAM 30-m telescope to map out 
emission from several spectral lines in order to under- 
stand the kinematic and spectral structure of TMC-1C. 
In November 1998, we mapped the spectral line maps 
of the C 17 O(l-0), C 17 0(2-l), C 18 0(2-l), C 34 S(2~1), 
DCO+(2-l), DCO+(3-2), N 2 H+(l-0) transitions. The 
inner 2' of TMC-1C were observed with 20" spacing in 
frequency-switching mode, and outside of this radius the 
data were collected with 40" sampling. The data were 
reduced using the CLASS package, with second-order 
polynomial baselines subtracted from the spectra. The 
system temperatures, velocity resolution, beam size and 
beam efficiencies are listed in Table [TJ 

3. RESULTS 
3.1. Spectra 

The spectra taken at the peak of the dust column den- 
sity map are shown in Figure [TJ The integrated intensity, 
velocity, line width and RMS noise for each transition is 
given in Table [2] From the figure it is evident that self- 
absorption is present everywhere, except in C 17 and 
C 18 lines. Clear signatures of inward motions (brighter 
blue peak; e.g. Myers et al. 1997, see also Sect. I4.3[) 
are only present in the high density tracers N2H + and 
DCO + , which typically probe the inner portion of dense 
cores (e.g. Caselli et al. 2002b; Lee et al. 2004). The 
C 34 S(2-1) line appears to be self-absorbed at the cloud 
velocity, but the spectrum is too noisy to confirm this. 

To highlight the extent of the "infall" asymmetry in 
N2H + (l-0), Fig. [2] shows the profile of the main hyper- 
fine component of the N2H + (l-0) transition (Fi,F = 2,3 
— » 1,2) across the whole mapped area (see Sect. 13. 2[) . 
We will discuss these spectra in more detail in Sect. 14.31 
but for now it is interesting to see how the profile shows 
complex structure, consistent with inward motion (red 
boxes) as well as outflow (blue boxes) and absorption 
from a static envelope. 

3.2. Maps 

From Fig.[T]and[2j it is clear that the main N 2 H + (l-0) 
hyperfine components are self-absorbed around the dust 
peak position. Therefore, a map of the N2H + (l-0) in- 
tensity integrated under the seven components will not 
reflect the N2H+ column density distribution. However, 
the weakest component (Fi F = 10 — > 11) is not af- 
fected by self-absorption, as shown in Fig. [3j where the 
weakest and the main (Fi F = 23 — > 12) components 
toward the dust peak position (the most affected by self- 
absorption) are plotted together for comparison. The two 
hyperfines in Fig.[3]have very different profiles: the main 
component is blue-shifted, suggestive of inward motions 
(see Sect. 14.3] ). whereas the weak component is symmet- 
ric and its velocity centroid is red shifted compared to 
the main component, indicating optically thin emission. 
Given that the self-absorption is more pronounced at the 
position of the dust peak, where the N2H + (l-0) optical 
depth is largest, we conclude that the weak component 
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is likely to be optically thin across the whole TMC-1C 
core. 

Thus, in the case of N2H + (l-0) self-absorption, we 
used the weak hyperfine component, divided by 1/27 
(its relative intensity compared to the sum, normal- 
ized to unity, of the seven hyperfines), to determine 
the N2H + (l-0) integrated intensity, line width, and, as 
shown in Sec. [SJ the N2H 4 " column density. In this anal- 
ysis, only spectra with signal to noise (S/N) for the weak 
component > 3 have been considered. In all other cases, 
N2H + (l-0) did not show signs of self-absorption and a 
normal integration below the 7 hyperfines has been per- 
formed. 

Based on hyperfine fits to the C 17 O(l-0) transition and 
a comparison of the relative strengths of the three com- 
ponents, we see that the C 17 O(l-0) emission is optically 
thin throughout TMC-1C. Although the noise is gener- 
ally too high in the C 17 0(2-l) data to make indisputable 
hyperfine fits, the results of such an attempt suggest that 
the C 17 0(2-l) lines are also optically thin, which is ex- 
pected for thin C 17 O(l-0) emission and temperatures of 
~10 K. Thus, the integrated intensity maps of C 17 lines 
will reflect the C 17 column density distribution. 

In order to estimate the optical depth of the C 18 0(2-l) 
lines, we compare the integrated intensity of C 18 0(2-l) 
to that of C 17 0(2-l). If both lines are thin, then the 
observed ratio should be equal to the cosmic abundance 
ratio, ^ib.i7 = [ 18 0]/[ 17 0] = 3.65 (| Wilson fc Roodlfl99l 
iPenziasI Il981| ). We observe that the ratio of the inte- 
grated intensities Ris,n = 2.3 ± 0.9, which corresponds 
to an optical depth of the C 18 0(2-l) line of ns — 1.5. 

Integrated intensity maps of C 17 O(l-0), C 17 0(2-l), 
C 18 0(2-l), DCO+(2-l) and N 2 H+(l-0) are shown in 
Fig. [U Note that the N 2 H + integrated intensity map 
peaks right around the position of the dust column peaks, 
which is not true for C 17 and C 18 0. We do not present 
integrated intensity maps of C 34 S(2-1) or DCO+(3-2), 
which have lower signal to noise. 

4. ANALYSIS. I. KINEMATICS 

4.1. Line Widths 

Ammonia observations have shown that TMC-1C is a 
coherent core, having a constant line width across the 
core at a value slightly higher than the thermal width, 
and increasing outside the " coherent" radius, ~0.1 pc 
(|Barranco fc Goodman! [1998D . Our N2H + observations 
of TMC-1C show that the line width remains constant, 
at a value ^2 times higher than the thermal line width, 
over the entire core (see Fig. [5] for a map of the line 
width and a plot of line width vs. radius), though the 
dispersion in the N2H + line width is very large. This re- 
sult is in agreement with the NH3 observations of TMC- 
1C, and is not consistent with the decreasing N 2 H + and 
N2D" 1 " linewidths at larger radii seen in L1544 and L694- 
2, which in other important ways (density and temper- 
ature structure, velocity asymmetry seen in N2H + (l-0)) 
closely resemble TMC-1C. To make sure that the lack 
of correlation of the N2H + (l-0) line width is not due 
to geometric effects, considering the elongated structure 
of TMC-1C, we also plotted the N 2 H+(l-0) line width 
as a function of antenna temperature, and found similar 
results. 

This behavior may be due to the coherence of the cen- 



tral portion of the core, which has nearly constant length 
along the line of sight, and thus the velocity dispersion 
comes from regions of the core that have similar scales 
(see Fig. 4 in Goodman et al. 1998). Cores formed 
by compressions in a supersonic turbulent flow naturally 
develop these region s of constant length at their centers 
(|Klessen et all I2005D . Another reason that the N 2 H+ 
line widths appear constant across the cloud could be 
the different "infall" velocity profile, with the velocity 
peaking farther away from the dust peak than in L1544 
and L694-2 (though the projected velocity would still 
be at its maximum at the dust peak). In the case of 
L1544. ICaselli et all (l2002aft show ed that the N2H+(l-0) 
line profile is consistent with the ICiolek fc Basul (2000) 
model at a certain time in the cloud evolution, where the 
"infall" velocity profile peaks at a radius of about 3,000 
AU (see also Myers 2005 for alternative models with sim- 
ilar radial velocities). Indeed, in Sec. 14.31 we show that 
the extent of the asymmetry seen in N2H + (l-0) suggests 
that the peak of the inward motions is at about 7,000 AU, 
so that one does not expect to see broader line widths 
within this radius. In fact, the binned data in Fig.[5]show 
a hint of a peak at about 50" (7000 AU at the distance 
of Taurus) . 

In order to compare the thermal and non-thermal line 
widths in TMC-1C, we assume that the gas temperature 
is equal to 10 K and use the formulae: 

y V ^ m H J 

Av NT = V(^obs) 2 - (A«t) 2 (2) 

where /x is the molecular weight of the species and uih 
is the mass of hydrogen. 

Figure [6] shows the non-thermal line width plotted 
against the thermal line width at the position of the dust 
column density peak. To allow a fair comparison, all 
the data in the figure have been first spatially smoothed 
at the same resolution of the NH3 map (1'). Although 
C 17 0, C 18 0, DCO+ and N 2 H+ all have similar molec- 
ular weights, they have significantly different values for 
their non-thermal line widths. The thermal line width 
is much smaller than the non-thermal line width for the 
molecules C 17 and C 18 0, while the ratio is closer to 
unity for DCO + and N2H" 1 ". This suggests that the iso- 
topologues of CO are tracing material at larger distances 
from the center, with a larger turbulent line width, than 
are DCO + and N 2 H + , which presumably are tracing the 
higher density material closer to the center of TMC-1C. 
The C 18 0(2-l) line is slightly thick, and this is probably 
the reason of its slightly larger line width when compared 
to the thin C 17 0(2-l) line, as shown in Figure [BJ For 
each transition observed, we see no clear correlation be- 
tween the observed line width and the thermal line width 
(and therefore with temperature, column density and dis- 
tance from the peak column density, see Section [5~2j) . As 
in the study of depletion, the lower signal to noise in 
C 34 S(2-1), DCO+(2-l) and DCO+(3-2) make any pos- 
sible trends between Av and Avthermai more difficult to 
determine. From Fig. [6] we note that NH3 and N2H + 
have similar non-thermal line widths, which makes sense 
given that N2H 4 " and NH3 are expected to trace similar 
material (e.g. Benson et al. 1998). However, this result 
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is in contrast with the findings of T afalla et al.l (|2004f ) 
who found narrower NH3 line widths towards L1498 and 
L1517B. 

We finally note that the line widths that we measure 
in C 17 and C 18 are larger than the N2H+ lincwidths 
throughout TMC-1C, which co ntrasts with the results 
seen in C 18 and N 2 H+ in B68 (|Lada et al.ll2003l) . This 
is consistent with the fact that TMC-1C, unlike B68, 
is embedded in a molecular cloud complex and it is not 
an isolated core. Thus, CO lines in TMC-1C also trace 
the (lower density and more extended) molecular mate- 
rial, part of the Taurus complex, where larger ranges of 
velocities are present along the line of sight. 

4.2. Velocity Gradients 

In order to study the velocity field of TMC-1C, we de- 
termine the centroid velocities for C 17 0(2-l), C 18 0(2- 
1), C 34 S(2-1), DCO+(2-l) and DCO+(3-2) with Gaus- 
sian fits. The centroid velocities of the C 17 0(1 — 0) and 
N2H + (l-0) lines are determined by hyperfine spectral 
fits. For those N2H + spectra that show evidence of self- 
absorption, the velocity is derived from a Gaussian fit to 
the thinnest component. The velocity gradient at each 
position is calculated by fitting the velocity field with the 
function: 

dv dv 
vi sr = v Q + —Aacosd + —ASsm6 (3) 
ds ds 

where v is the bulk motion along the line of sight, Aa 
and AS are RA and DEC offsets from the position of 
the central pixel, dv/ds is the magnitude of the velocity 
gradient in the plane of the sky, and 9 is direction of 
the velocity gradient. The fit to the velocity gradient is 
based on fitting a plane through the posi tion-position ve- 
locity cube as in I Goodman et alj (119931) ffor the "total" 
gradient across the cloud) and in iCaselli eTaLl (|2002bf ) 
(for the "local" gradient at each position). The fit for 
the "total" velocity gradient gives a single direction and 
magnitude for the entire velocity field analyzed. The 
"local" velocity gradient is calculated at each position in 
the spectral line maps based on the centroid velocities of 
the center position and its nearest neighbors, with the 
weight given to the neighbors decreasing exponentially 
with their distance from the central position. 

Analysis of ammonia observations with ~ 60" reso- 
lution of TMC-1C indicate an overall velocity gradi- 
ent of 0.3 km s -1 pc~^ dire cted 129 degrees East of 
North (|Goodman et al.lll993D . The velocity field that 
we measure in TMC-1C has spatial resolution three 
times greater (~20") than the ammonia study, and re- 
veals a pattern more complicated than that of solid body 
or differential rotation. The velocity fields measured 
by C 17 O(l-0), C 17 0(2-l) and C 18 0(2-l) are shown in 
Fig. [7] Although there is a region that closely resembles 
the velocity field expected from rotation (gradient arrows 
of approximately equal length pointing in the same di- 
rection), the measured velocities vary from blue to red 
to blue along a NW - SE axis. The N 2 H+(l-0) veloc- 
ity fields (shown in Fig. O also follow the same blue to 
red to blue pattern along the NW - SE axis, but the ob- 
servations cover a somewhat different area than the CO 
observations, which complicates making a direct compar- 
ison. Taken as a whole, it is clear that there is an ordered 



velocity field in portions of the TMC-1C core, and that 
the lower density CO tracers "see" a velocity field similar 
to that probed by N 2 H + lines, which trace higher density 
material. In any ca se, the velocity field th at looks like 
rotation reported in iGoodman et al.l (|1993h turns out to 
be more complicated when seen over a larger area with 
finer resolution. The direction and magnitude of the ve- 
locity gradient in the region that resembles solid body 
rotation is shown in Fig. [8] for each transition. 

4.3. Inward Motions 

To quantify the velocity of the inward motion from 
the N 2 H+(l-0) line across the TMC-1C cloud, we use 
a simple two-layer model, similar to that described by 
iMvers et al.1 (Jl996) . This model assumes that the cloud 
can be divided in two parts with uniform excitation tem- 
perature (T cx , gradients in T ex between the two layers as 
in De Vries & Myers 2005, are not considered here), line 
width (Av), optical depth (r) and LSR velocity (Vlsr) 
and that the foreground layer has a lower excitation tem- 
perature. For simplicity, we also assume that the seven 
hyperfines have the same T cx , which is a very rough as- 
sumption in_j^gionsof large optical depth, as recently 
found bv I Daniel et al.l (|2006h . Despite of the simplicity 
of the model, we find good fits to the seven hyperfine 
lines and determine the value of the velocity difference 
between the two layers, which can be related to the "in- 
fall" velocity. 

In Fig. [S] we present five spectra which represent a cut 
across the major axis of the core, passing through the 
dust peak. For display purposes, the spectra have been 
centered to velocity, subtracting the LSR velocity ob- 
tained from a Gaussian fit to the weak hyperfine compo- 
nent (for offsets [-40,60], [-20,40], [0,20], where the self- 
absorption in present) or hfs fits in CLASS 7 (for offsets 
[20,0] and [40,-20]). The Vlsr velocity is shown in the top 
right of each panel. The cut is from South-East (offset 
[40,-20], see Fig. [2]) to North- West (offset [-40,60]). The 
first thing to note in the figure is that clear signs of self- 
absorption and asymmetry are present toward the dust 
peak and in the North- West, but not in the two Southern 
positions. This trend can also be seen as a general fea- 
ture in Fig. [21 where it is evident that asymmetric lines 
are more numerous North- West of the dust peak. 

The excitation temperature, total optical depth, line 
width and the velocity (V — Vlsr, see Fig. |TJ]) of the fore- 
ground (F) and background (B) layers are reported in Ta- 
ble^ To find the best fit parameters, we first performed 
an hfs fit to the [20,0] spectrum, which is the closest spec- 
trum to the dust peak not showing self-absorption. The 
values of ttot, T cx and line width obtained from this fit 
have been adopted for the background emission at the 
dust peak position and the best fit has been found by 
adding the foreground layer and minimising the residu- 
als. For the two spectra North-West of the [0,20] po- 
sition, adjustment to the parameters of the background 
layer were necessary to obtain a good fit. We point out 
that the five spectra we chose for this analysis are rep- 
resentative of the whole area surrounding the TMC-1C 
dust peak, where a mixture of symmetric, blue-shifted 

7 CLASS is part of the Grenoble Image and 
Line Data Analysis Software (GILDAS), available at 
|http://www.iram .fr/IRAMFR/GILDAS/gildas.htmi). 
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and red-shifted spectra are present. As already stated, 
the majority of the asymmetric spectra show inward mo- 
tions and extend over a region with radius ^7000 AU 
(see Fig.©. 

The properties derived for the foreground layer (T ox 
~3.3-3.5 K, r TO T ^10-15, and Av ~0.2 kms' 1 ) have 
been used as input parameters in a Large Velocity Gra- 
dient (LVG) code 8 for a uniform medium and found to be 
consistent with the N2H + (l-0) tracing gas at a density 
n(H2) — 5xl0 3 cm -3 , kinetic temperature Tki n — 10 K 
and with column density vV(N2H + ) ~ 5xl0 12 cm -2 , val- 
ues comparable to those found for the background layer 
(see Table E and Sec. E3) . 

It is interesting that the maximum of the line-of- 
sight component of the inward velocity (~0.15 kms - ) 
is found toward the dust peak, whereas one pixel away 
from it, the inward velocity drops to 0.05 kms' 1 . This 
is suggestive of a geometric effect, in which the inward 
velocity vector is directed toward the dust peak, so that 
only a fraction cos(9) (with 8 the angle between the l.o.s. 
and the infall velocity direction) of the total velocity is 
directed along the line of sight in those positions away 
from the dust peak. Of course, our simplistic model pre- 
vents us to go further than this, i.e. the uncertainties 
are too large to build a 3D model of the velocity profile 
within the cloud. 

As shown in Fig. [21 in the North-West end of the 
TMC-1C core (around offset [-250,150]), there are other 
signatures of inward motions, which may indicate the 
presence of another gravitational potenti al well. This 
sugges tion is indeed reinforced by Fig. 8 of ISchnee et al.1 
(2007), which shows high extinction and low tempera- 
tures in the same direction. Unfortunately, the contin- 
uum coverage is not good enough to attempt a detailed 
analysis, but it appears evident that the extension to- 
ward the North- West is another dense core connected to 
the main TMC-1C condensation with lower density and 
warmer gas and dust. 

5. ANALYSIS. II. COLUMN DENSITY AND TEMPERATURE 

5.1. Gas Column Density 

To derive the column density of gas from each 
molecule, we assume that all rotation levels are char- 
acterized by the same excitation temperature T ex (the 
CTEX method, described in Caselli et al. 2002b). In 
case of optically thin emission, 

8irW gi 1 1 



and 



Ej = J(J+l)hB 



(6) 



and B is the rotational constant (see Table 0] for the 
values of the constants). 

The C 17 O(l-0) and (2-1) lines have hyperfme struc- 
ture, enabling the measurement of the optical depth. We 
find that the lines are optically thin throughout the core. 
To determine the column density, we assume an excita- 
tion temperature of 11 K, which is the average value of 
T cx found from our C 18 data around the dust peak po- 
sition, as explained in Sect. 15.31 In the case of C 18 0(2-l) 
lines, we correct for optical depth before determining the 
column density, using the correction factor: 



1-e- 



(7) 



Ni 



As explained in Sect. I3.1] N2H + (l-0) lines show clear 
signs of self-absorption in an extended area around the 
dust peak. To determine the column density across the 
core, first we select spectra without self-absorption, and 
those with high S/N (i.e. with W/a w > 20, with W = 
integrated intensity; see Caselli et al. 2002b) have been 
fitted in CLASS to find T ex and r. The mean value of T EX 
found with this analysis (4.4±0.1 K) has been used for 
all other positions where an independent estimate of T ox 
was not possible (i.e. for self-absorbed or thin lines). For 
optically thin N2H + (l-0) transitions, the intensity was 
integrated below the seven hyperfme and the expression 
([4]) used to determine the total column density. 

In cases of self-absorbed spectra, the N2H + column 
density has been estimated from the integrated intensity 
of the weakest (and lowest frequency) hyperfme compo- 
nent (FiF = 1 0— >1 1), using eq. [4] (assuming T ex =4.4 K) 
and multiplying by 27 (the inverse of the hyperfme rel- 
ative intensity). The weakest component is not affected 
by self-absorption, as shown in Fig.[3l suggesting that its 
optical depth is low. We have checked that these two dif- 
ferent methods approximately give the same results by 
measuring the N 2 H + column density with both proce- 
dures in those cases where self-absorption is not present 
and where the weakest hyperfme component has a S/N 
ratio of at least 4. We found that the two column density 
values agree to within 10%. 

DCO + (2-l) and DCO + (3-2) lines are clearly self- 
absorbed and the column density determination is very 
uncertain (given that there are no clues about their op- 
tical depth and excitation temperature). The estimates 
Qstod in Table [5] should be considered lower limits. T ox 



X 3 A g u J v {T ex ) - J v {T bg ) 1 - exp(-/w/fcT e:c ) gi exp 

(4) 

where A and v are the wavelength and frequency of the 
transition, k is the Boltzmann constant, h is the Planck 
constant, A is the Einstein coefficient, gi and g u are the 
statistical weights of the lower and upper levels, J v (T ex ) 
and J u (Tb g ) are the equivalent Ray leigh- Jeans excitation 
and background temperatures, W is the integrated in- 
tensity of the line. The partition function (Q ro t) and 
the energy of the lower level (Ei) for linear molecules are 
given by: 



Qrot = ^ (2 J + !) eM-Ej/kT) 



(5) 



.7=0 



' available at http:/ /www. strw.leidenuniv.nl/~moldata/radex.php 



(^Ei/4c% x has beed assumed, based on the fact that the 
DCO+ lines are expected to trace similar conditions than 
N2H + . C 34 S(2-1) spectra have low sensitivity, and the 
lines are affected by self-absorption (see Fig. QJ, so the 
derived C S column density is highly uncertain. 

The N 2 H+ abundance, A(N 2 H+) (= 
iV(N2H+)/A(H 2 )), toward the dust peak is 1.6xl0- 10 , 
identical (within the errors) to that de rived toward 
the L1544 dust peak (ICrapsi et all l2005l e.g.). This 
is interesting considering that tV(H2) in TMC-1C is 
1.6 times lower than in L1544, in which N2H+ closely 
follows the dust column (as already found in previ- 
ous work). However, unlike L1544, where the N 2 H + 
abundance appears constant with impact parameters 
(e.g. Tafalla et al. 2002 and Vastel et al. 2006), in 
TMC-1C the N 2 H + abundance increases away from 



() 
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the dust peak by a factor of about two within 50 , as C O line, which is coincident with the kinetic temper 



shown in Fig. [TO] (see also Fig. [H in Sect. ETIj) . Fig. [TO] 
displays Ay (see Sect. HT2J), A(N 2 H+) and A(N 2 H+), 
normalized to the corresponding maximum values (63.2 
mag, l.lxlO 13 cm" 2 , and 4.8xl0~ 10 , respectively) in 
two cuts (one in right ascension and one in declination) 
passing through the dust peak. One point to note is 
that the abundance derived at the dust peak (marked by 
the black dotted line) is the minimum value observed, 
indicating moderate (factor of ~2) depletion. 

5.2. Dust Column Density 

In lSchnee et "all (pool we used SCUBA and MAMBO 
maps at 450, 850 and 1200 ^m to create column den- 
sity and dust temperature maps of TMC-1C. In this pa- 
per we smooth the dust continuum emission maps to the 
20" spacing of the IRAM maps and then derive Ay and 
Td to facilitate a direct comparison of the gas and dust 
properties. At each position, we make a non-linear least 
squares fit for the dust temperature and column density 
such that the difference between the predicted and ob- 
served 450, 850 and 1200 /im observations is minimized. 
The errors a ssociated with such a fitting procedure are 
described in ISchnee et al.1 (j2007) . Dust column density 
a nd temperature map s of TMC-1C are shown in Fig. 8 
of ISchnee et all (120071). 

From Fig. 8 in ISchnee et"afl (|2007l ). it is clear that 
there is an anti-correlation between extinction and dust 
temperature (as also predicted by theory, e.g. Evans et 
al. 2001, Zucconi et al. 2001, Galli et al. 2002). To bet- 
ter show this, the two quantities are plotted in Fig. [TT1 
The data in Fig. [IT] are not smoothed to the IRAM 30m 
beam at 3 mm, since this is only a dust property in- 
tercomparison and does not refer to the gas properties. 
Higher Ay (60 < Ay < 90 mag) and lower dust temper- 
atures (5 < Idust < 6 K) are detectable at the higher 
resolution. We compare the Td — Ay relationship seen 
in TMC-1C with that predicted bv lZucconi et all (|2001l ) 
for an externally heated pre-protostellar core (the solid 
red line in Fig. fTTj). We find that at high column density 
(Ay > 30) the observed dust temperature in TMC-1C is 
lower than that of the model core, while at low column 
density (10 < Ay < 20) the observed dust temperature 
is higher than the model predicts. However, given that 
the model predicts the dust temperature at the center of 
a spherical cloud, and that the geometry of TMC-1C is 
certainly not spherical, only a rough agreement between 
the model and observations should be expected. 

5.3. Gas Temperature 

Because of its low dipole moment, CO is a good gas 
thermometer, given that it is easily thermalized at typ- 
ical core densities. However, it is now well established 
that CO is significantly frozen onto dust grains at densi- 
ties > 10 5 cm -3 (one exception being L1521E; Tafalla & 
Santiago, 2004b) and this is also the case in TMC-1C. 
Therefore, at the dust peak we do not expect to mea- 
sure a gas temperature from CO of ~7 K, but instead 
a higher value reflecting the temperature in the outer 
layers of the cloud. The lines available for this anal- 
ysis are: C 17 O(l-0), C 17 0(2-l), and C 18 0(2-l). The 
C 17 0(2-l)/C 18 0(2-l) brightness temperature ratio has 
been used to derive the excitation temperature of the 



ature if the line is thermalized. To test the hyposthe- 
sis of thermalization we use an LVG (Large- Velocity- 
Gradient) program to determine at which volume den- 
sity and kinetic temperature the observed C 17 O(l-0) 
and C 17 0(2-l) brightness temperatures can be repro- 
duced. We use a one-dimensional non -LTE radia- 
tive transfer code (Ivan der Tak et~aTll2007ft available at 
|http: / / www.strw.leidenuniv.nl~ moldata/radex.html] 

5.3.1. C 17 0(2-1) and C 18 0(2-1) as a measure o/T kin 

These two lines have similar frequencies, so the corre- 
sponding angular resolution is almost identical and no 
convolution is needed. Following a similar analysis done 
with the J =1-0 transition of the two CO isotopologues 
(Myers et al. 1983), the optical depth of the C 18 0(2-l) 
line (rig) can be found from: 



T mb [C 18 Q(2-l)] _ i-e-^s 
T mb [C 17 O(2-l)]- d - b0X tib : 



(8) 



where T mb [i] is the main beam brightness temperature 
of transition i (assuming a unity filling factor) . The last 
term in the right hand side is the optical depth correc- 
tion which is used to determine the total column density 
of C 18 in a plane parallel geometry, which most likely 
applies to CO emitting regions, i.e. the external core 
layers. 

Once rig is measured, the excitation temperature (T ex ) 
of the corresponding transition (thus the gas kinetic tem- 
perature, if the line is in local thermodynamic equilib- 
rium) can be estimated from the radiative transfer equa- 
tion: 



T mb - [J„(T«) - J„(T bg )](l - e- r ), 



(9) 



where J v {T ex ) and J v (Tb g ) are the equivalent Reyleigh- 
Jeans temperatures, with 



J AT) = 



exp(T /T) - 1 : 



(10) 



To = hv/k&, and v the frequency of the C 18 0(2-l) line 
(see Table P. 

Figure [12] (left panel) shows the results of this analysis. 
The set of data points in Fig. [T^] is limited to only those 
spectra with T mb [C 17 0(2--l)]/T rms > 10, to avoid scat- 
ter due to noise. The error associated with the gas tem- 
perature has been calculated by propagating the errors 
on rig and T mb [C ls O(2-l)] into Eq.[9]and its expression 
is given in the appendix. 

It is interesting to note that the C 18 0(2-l) excitation 
temperature appears to decrease away from the center, 
and, if this line is thermalized, it suggests that the kinetic 
temperature also drops away from the center, in contrast 
with the dust temperature. Indeed, the two quantities 
are completely uncorrelated in the 12 common positions 
with large S/N C 18 0(2-l) spectra (not shown). T GX is 
close to 12 K at the core center, whereas it drops to 9- 
10 K one arcmin away from the dust peak. Is this drop 
du e to a decreas i ng ga s temperature, as recently found 
by iBergin et ail (|2006D in the Bok Globule B68? Un- 
like B68, we believe that our result is due to the volume 
density decrease. Indeed, the critical densities of the J 
= 2-1 lines of C 17 and C ls O are a few x 10 4 cm" 3 , 
so that only if the volume density traced by one of the 
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two isotopologues is larger than, say, 5x10 cm -3 can 
the J = 2-1 lines be considered good gas thermometers. 
In the next subsection, we investigate this point more 
quantitatively. 

5.3.2. C 17 O(l~0) and C 17 0(2-1) to measure 
T ox /e 17 0(2-1)] 

In Fig. Q2] (right panel), the brightness temperature 
ratio of the C 17 O(l-0) and C 17 0(2-l) lines is plotted 
as a function of distance from the dust peak. Because 
of the different angular resolutions at the 2-1 and 1-0 
frequencies, the 1 mm data have been smoothed to the 
3 mm resolution and both data cubes have then been 
regridded, to allow a proper comparison. The ratio is in- 
deed increasing towards the edge of the cloud, consistent 
with our previous finding of a T ox [C 18 0(2-l)] drop in the 
same direction (see left panel). 

Both C 17 O(l-0) and (2-1) lines possess hyperfmc 
structure, which provides a direct estimate of the line 
optical depth. Using the hfs fit procedure available in 
CLASS, we found that all over the TMC-1C cloud, both 
lines are optically thin. This means that it is not possible 
to derive the excitation temperature in an analytic way, 
so we use the LVG code introduced in Sect. 15.31 This 
code assumes homogeneous conditions, which is likely to 
be a good approximation for the region traced by CO iso- 
topologues. In fact, because of freeze-out, CO does not 
trace the regions with densities larger than about 10 5 
cm -3 (see below and Sect. I6.1|l . so that the physical con- 
ditions traced by CO around the dust peak are likely to 
be close to uniform (n(H2) ~ a few times 10 4 and about 
constant temperature) . This is also supported by the in- 
tegrated intensity CO maps, which appear extended and 
uniform around the dust peak (see Fig. [4]). 

To better understand this result, the LVG code has 
been run to see how changes in volume density and 
gas temperature affect the line ratio. This is shown in 
Fig. [13l where the top panel shows the T m b ratio as 
a function of T^ n for a fixed value of the volume den- 
sity (n(H 2 ) = 2xl0 4 cm -3 ), a C 17 Q column density of 
10 15 cm -2 (as found in Section [BTTjl . and a line width of 
0.4km s - , as observed. The horizontal dashed lines en- 
close the range of T m b ratios observed in TMC-1C and 
reported in Fig. rTJJ Thus, the observed T m b range (at 
this volume density) corresponds to a range of gas tem- 
perature between 11 K toward the dust peak and ~7 K 
away from it, thus confirming our previous findings of a 
decreasing C 18 0(2-l) excitation temperature away from 
the dust peak. 

In the bottom panel of Fig. [131 the same brightness 
temperature ratio is plotted as a function of n(H2), for 
a fixed kinetic temperature (Tki n =H K), iV(C 17 0) = 
10 15 cm -2 and Aw=0.4 kms - , as before. The black 
curve shows this variation and, not surprisingly, the ob- 
served range of T m b ratios can also be explained if the 
volume density (traced by the C 17 lines) decreases from 
~4x 10 4 cm -3 toward the dust peak to ~ 6 x 10 3 cm -3 
away from it (the point farthest away being at a pro- 
jected distance of 80", or 11,000 AU, see Fig. [12). Note 
that the volume density traced by the C 17 line towards 
the dust peak is significantly lower than the central den- 
sity of TMC-1C (~5xl0 5 cm" 3 ; see Schnee et al. 2007), 
once again demonstrating that CO is not a good tracer 



of dense cores. The bottom panel of Fig. [12] is consis- 
tent with a volume density decrease away from the dust 
peak, or, more precisely, a lower fraction of (relatively) 
dense gas intercepted by the C 17 lines along the line 
of sights. In the same plot, the red curves show the ex- 
citation temperatures of the C 17 O(l-0) and (2-1) lines 
vs. n(H 2 ). Note that T ox [C 17 0(2-l)] = T kin only when 
the density becomes larger than ~10 5 cm" 3 . Thus, the 
C 17 0(2-l) line is sub-thermally excited in TMC-1C. 

In summary, the rise in the C 17 O(l-0)/(2-l) brigth- 
ness temperature ratio away from the dust peak ratio 
can be caused by either a gas temperature decrease or a 
volume density decrease (or both) . Considering that the 
dust (and likely the gas; see the recent paper by Crapsi 
et al. 2007) temperature is clearly increasing away from 
the dust peak, we believe that the drop in T ex observed 
both using the C 18 and C 17 lines is more likely due 
to a drop in the volume density traced by these species. 
This is reasonable in the case of a core embedded in a 
molecular cloud complex, such as TMC-1C, where the 
fraction of low density material intercepted along the line 
of sight by C 18 and C 17 observations is significantly 
larger than in isolated Bok Globules such as B68 (see 
Bergin et al. 2006). In any case, a detailed study of 
the volume density structure of the outer layers of dense 
cores will definitely help in assessing this point. 

6. ANALYSIS. III. CHEMICAL PROCESSES 

6.1. Molecular Depletion 

By comparing the integrated intensity maps of CO 
isotopologues and N 2 H+(l-0) (Fig. [4]) in TMC-1C with 
the column density implied by dust emission (Fig. 8 in 
iSchnee et all (f2007h L we sec that at the location of the 
dust column density peak the CO emission is not peaked 
at all. The N2H + (l-0) emission peaks in a ridge around 
the dust column density maximum, not at the peak, but 
in general N2H 4 " traces the dust better than the C 18 
emission does. Below, we measure the depletion of each 
observed molecule and compare our results to similar 
cores. 

Previous molecular line observations of starless cores 
such as L1512, L1544, L1498 and L1517B, consistently 
show that C O and its isotopologues are significa ntly de- 
pleted, e.g. (|Lee et all [20031 : iTafalla et al.ll2004D . How- 
ever, other molecules such as DCO + and N 2 H + are typ- 
ically found to trace the dust emission well, (e.g. Caselli 
et al. 2002b; Tafalla et al. 2002, 2004), although there is 
some evidence of their depletion in the center of chem - 

2002) . 

2003) . 



ically evolved cores, such as B68 (IBergin et al 
L1544 (jCaselli et al.l l2002b) and L1512 (|Lee et al.l 



In order to measure the depiction in TMC-1C, we define 
the depletion factor of species i: 



f D (i) = X can (i) 



N(H- 



2 J dust 



N(i) 



(ii) 



where X can is the "canonical" (or undepleted) fraction 
abundance of species i with respect to H2 (see Table [4]) , 
N{H2)dust is the column density of molecular hydrogen 
as derived from dust emission, and N(i) is the column 
density of the molecular species as derived in Section [STTl 
The derived depletion factors for each molecule (except 
for DCO + and C 34 S, where the column density determi- 
nation is quite uncertain, as explained in Sec. 15. ip in 
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each position with signal to noise > 3 are plotted against 
the dust-derived column density in Figure [TU The typ- 
ical random error in the derived depletion is shown in 
each panel, and is derived from the noise in the spec- 
tra. Uncertainty in the derived column density from 
dust emission is dominated by calibration uncertainties 
in the bolometer maps, and is not included in this cal- 
culation, nor is the uncertainty in the calibration of the 
spectra (~20%), which would adjust the derived deple- 
tion factors systematically, but would not alter the ob- 
served trends. Depletions factors are found to increase 
with higher dust column density (Fig. [T4J) . In TMC-1C 
we see a linear relationship between C 17 and C 18 de- 
pletion and dust-derive d column den s ity, w hich has also 
been seen in C 18 by ICrapsi et all ((2001 in the core 
L1521F, which cont ains a Very Low Luminosity Object 
(|Bourke et al.l 120061 ). To check the impact of resolu- 
tion on the derived depletion, we compare the C 18 0(2- 
1) depletion when smoothed to 20" (the resolution of 
the C 17 O(l-0) data) with that derived from smoothing 
C 18 0(2-l) to 14" (the resolution of the bolometer data). 
We find no systematic difference between the two calcu- 
lations of the depletion, and a 13% standard deviation in 
the ratio of the derived depletions. 

The depletion factor and column density at the posi- 
tion of the dust peak is listed in Table [5] for each tracer. 
The depletion factor measured in N2H + clearly follows a 
different trend compared to the CO isotopologues. First 
of all, in Fig. [14] the N2H 4 " depletion factor is allowed to 
have values below 1, because of our (arbitrary) choice for 
the "canonical" abundance of N2H + assumed here to be 
equal to 1.4xl0~ 10 , the average value across TMC-1C. 
We point out that a "canonical" abundance for N2H 4 " is 
much harder to derive than for CO, because N2H 4 " lines 
are much harder to excite (and thus detect) in low den- 
sity regions where depletion is negligible. Nevertheless, 
Fig. 20 shows that the N2H" 1 " depletion factor monotoni- 
cally increases (as in the case of CO) for N(H2) > 3 x 10 22 
cm~ 2 (or Ay > 30 mag). This is clear evidence o/N2H + 
depletion in the core nuclei, in a central region with ra- 
dius ^6000 AU, where A\j > 30 mag (see also Fig. [TO]) . 

The dispersion in the N2H + depletion factor vs. N(H2) 
relation is very large with no obvious trend at lower 
Ay values (N(H 2 ) < 3 x 10 22 cm -2 ), which we believe 
is due to our choice of excitation temperature where 
the N2H + (l-0) line is optically thin. As explained in 
Sect. 15. ll in the case of optically thin lines, T ox has been 
assumed equal to 4.4 K, the mean value derived from the 
optically thick spectra which do not show self-absorption 
(and which trace regions with Ay > 20 mag). There- 
fore, in all positions with N(H2) below 2xl0 22 cm -2 , 
where the volume density is also likely to be low, the 
assumed N2H + (l-0) excitation temperature is likely to 
be an overestimate of the real T cx . To see if this can 
indeed be the cause of the observed scatter, consider a 
cloud with kinetic temperature of 10 K, volume density 
of 3xl0 4 cm -3 , line width of 0.3 km s _1 and an excita- 
tion temperature of 4.4 K for the N2H + (l-0) line. Using 
the RADEX LVG program, this corresponds to a N2H+ 
column density of 10 12 cm~ 2 . If the density drops by 
a factor of two (whereas all the other parameters are 
fixed), the N2H + (l-0) excitation temperature drops to 
3.6 K. In these conditions, using T cx = 4.4 K instead of 
3.6 K, in our analytic column density determination (see 



Sect. 15. ip . implies underestimating N(N2H + ) by 50%. 
Therefore, our assumption of constant T ex can be the 
main cause of the observed fn scatter at low extinctions. 

Because of the anti-correlation between dust temper- 
ature and column density (see Fig. [TTj) . we expect that 
there will also be an anti-correlation between the deple- 
tion factor, frj, and dust temperature. Figure [T4l (right 
panels) shows the depletion factor for each molecule plot- 
ted against the line-of-sight averaged dust temperature. 
As expected, the depletion is highest in the low temper- 
ature regions though the lower signal to noise in C 34 S 
and DCO + make this somewhat harder to see. The anti- 
correlation between the depleti on factor and dust te m- 
perature has also been seen by iKramer et al.l (|1999f ) in 
IC5146 in C 18 0, though in TMC-1C the temperatures 
are somewhat lower. 

Our data clearly suggest that there is an increasing 
depletion of N2H + with increasing H2 column (and vol- 
ume) density. In previous work (e.g. Tafalla et al. 2002, 
2004; Vastel et al. 2006), the observed N 2 H+ abundance 
appears constant across the core, although the data are 
also consistent with chemical models in whi ch the N2H" 1 " 
abund a nce decreases by fa ctors of a few (jCaselli et al.1 
l2002bf) . iBergin et all (|2002f) also deduce small depletion 
factors for N 2 H + w hen comparing data to models and 
iPagani et al.l (|2005l ) found clear signs of N2H + deple- 
tions at densities above ^10 5 cm~ 3 . There is also ev- 
idence of N 2 H + dep letion towards the Class protostar 
IRAM 04191+1522 (jBelloche fc Andrei l200l in Taurus. 
The average N2H + abundance that we find in TMC-1C, 
relative to H 2 , is 1.4xl0~ 10 . 

What appears to be different from previous work is 
that the CO depletion factor towards the dust peak 
is relatively low (compared to, e.g., L1544), and, at 
the same time, N2H + (l-0) lines are bright over an 
extended region. If TMC-1C were chemically young 
(such as L1521E; Tafalla & Santiago 2004, Hirota et 
al. 2002), then there would be negligible CO freeze- 
out and low abundances of N2H + , given that N2H+ is a 
"late-type" molecule (i.e. its formation requires signifi- 
cantly longer times (factors >10) than CO and other C- 
bearing species). In TMC-1C we observe moderate CO 
and N2H + depletions, as well as extended N2H" 1 " emission 
with derived fractional abundances around 10~ 10 . To 
derive an approximate value of the average gas number 
density of the region where N2H + emission is present, we 
first sum all the observed N2H + (l-0) spectra (over the 
whole mapped area, with size ~450"xl70", correspond- 
ing to a linear geometric mean of about 40,000 AU; see 
Fig. [7]) and then perform an hfs fit in CLASS to derive the 
excitation temperature. We find T cx = 3.6±0.02 K and 
7V(N 2 H+) = 4.84±0.03 xlO 12 cm" 2 , which can be repro- 
duced with the LVG code if n(R 2 ) = 5xl0 3 cm~ 3 and 
Tkin = 11 K, as found in previous sections (see Sect. 14.31 
and I5.3|) . The average value of the extinction across the 
whole TMC-1C core is 23 mag, so that the corresponding 
N2H" 1 " abundance is 2xl0 -10 , close to the average value 
found before. How long does it take to form N2H + with 
fractional abundances of ~ 10~ 10 in region s with volume 
densities ~ 5 x 10 3 cm -3 ? iRoberts et all (|2004f) derive 
times > 3 x 10 5 yr at n(H2) ~ 10 4 cm~ 3 , so that this can 
be considered a lower limit to the age of the TMC-1C 
core. 

In summary, all the above observational evidence sug- 
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gests that the majority of the gas observed towards 
TMC-1C has been at densities ~10 4 cm -3 for at least a 
few times 10 5 yr and that material is accreting toward 
the region marked by the mm dust emission peak. We fi- 
nally note that the density profile of the region centered 
at the dust peak position is steeper (consistent with a 
power law; see Fig. 13 of Schnee & Goodman 2005) than 
found in other cores, so that Bonnor-Ebert spheres may 
not be the unique structure of dense cores in their early 
stages of evolution. 

6.2. Chemical Model 

Although TMC-1C is more massive than L1544 by 
a factor of about two, the physical structures of the 
two cores are similar: the central density of TMC-1C 
is ~ 5 x 10 5 cm -3 (factor of ~ 2 lower than L1544, 
according to Tafalla et al. 2002) and the central tem- 
perature is ^7 K, similar to the dust temperature de- 
duced by Evans et al. (2001) and Zucconi et al. (2001) 
in the center of L1544, and close to the gas tempera- 
ture recently measured by Crapsi et al. (2007), again 
toward the L1544 center. However, the chemical char- 
acteristics of the two cores appear quite different. In 
TMC-1C: 1. the observed CO depletion factor is about 
4.5 times smaller than in LI 544 (see Section 16.11 and 
ICrapsi et al.1 (|2005l ) ); 2. the deuterium fractionation is 
three times lower (|Crapsi et al.l 120051 ) than in L1544; 
3. the N 2 H + column density at the dust peak is two 
times lower and the C 17 Q column density is 1.7 times 
larger than in L1544 (jCaselli et all l2002bi ). All this is 
consistent with a younger chemical (and dyn amical) age 
(jShematovich et aUfeOOa lAikawa et al.ll2005h . 

To understand this chemical differentiation in objects 
in apparently similar dynamical phases, we used the sim- 
ple che mical model originally described in ICaselli et al. 
( 2002bh and more recently updated by IVastel et al. 
(2006). The model consists of a spherical cloud with 
density and tempe rature gradients as determined by 
ISchnee et all (|2007h . The model starts with H 2 , N 2 , CO 
and O in the gas phase, a gas-to-dust mass ratio of 100, 
and a Mathis, Rumpl & Nordsiek (1977; MRN) grain size 
distribution. Molecules and atoms are allowed to freeze- 
out onto dust grains an d deso rb via cosmic-ray impulsiv e 
heating (jCaselli et al.ll2002U lHasegawa fc Herbstlll993f ). 
The adopted binding energies of CO and N 2 are 1100 K 
and 982.3 K, respectively. The CO binding energy is in- 
termediate between the one measured for CO onto (i) icy 
mantles (1180 K; Collings et al. 2003 and Fraser et al. 
2004) and (ii) CO mantles (885 K; Oberg et al. 2005). 
The adopted value (1100 K) is the weighted mean of the 
two measured values, assuming that water is about four 
times more abundant than CO in the Taurus molecu- 
lar cloud (see Table 2 of Ehrenfreund & Charnley 2000 
and references therein) . See lOberg et al.l (|2005h for ad- 
sorption onto icy mantles. For the atomic oxygen bind - 
ing energy we used 750 K, as in IVastel et alJ (2006). 
The following para meters have also been assumed from 
IVastel et al.1 (|2006): (i) the cosmic ray ionization rate 
(1.3xl0~ 17 s" 1 ); (ii) the minimum size of dust grains 
(flmin = 5xl0 -6 cm ); (iii) the "canonical " abundance of 
CO (9.5 xl0~ 5 , fromlFrerking et al.l (I1982T): Civ) the stick- 
ing co efficient (5=1, as recently found by Bisscho p et all 
(2006) for CO and N 2 ); (v) the initial abundance of 



N 2 equal to 4xl0 -5 , i.e. about 50% the total abun- 
dance of nitrogen observed in the interstellar medium 
(|Mever et al.| [l997): (vi) the initial abundance of "met- 
als" (M+, in Fig. [I5J of 10~ 6 (fromlMcKeJ lfl989h ): and 
(vii) the initial abundance of Oxygen, fixed at a half the 
canonical abundance of CO (i.e. 13 times lower than the 
cosmic abundance; Meyer et al. 1998). 

The model is run until the C 17 column density to- 
ward the center of the cloud reaches the observed value 
(t = 8xl0 3 yr). During this time, the abundance of 
molecular ions is calculated within the cloud using steady 
state chemical equations with the instantaneous abun- 
dances of the ne utral species. To determine x(e ), the re- 
action scheme of lUmebavashi fc Nakanol (|1990l ) is used, 
where the abundance of the generic molecular ion "mH + " 
(essentially the sum of HCO+, N 2 H +, H 3 Q+ and their 
deuterated forms) is calculated (see lCaselli et al.1 J2002b) 
for more details). The calculated abundance profiles of 
the various species have then been convolved with the 
HPBWs of the 30m antenna at the corresponding fre- 
quencies and the derived column densities are in very 
good agreement with the observed quantities (within fac- 
tors of 2 for N 2 H + , N 2 D + and, of course, CO isotopo- 
logues), which is very encouraging, considering the sim- 
plicity of the model. 

The best-fit chemical structure of TMC-1C, reached 
after 10,000 yr, is shown in the left panel of Fig. [T5J 
Note that despite the similar binding energies of CO and 
N 2 , the HCO and DCO + drops are steeper than those 
of N 2 H + and N 2 D + , which is due to the fact that the 
CO freeze-out (although lower than in L1544) enhances 
the N 2 H + produ ction rate, as pointe d out by previous 
chemical models (Aik awa et al.l 12001) . Finally, we note 
that the CO depletion factor within the cloud, -Fb(CO) 9 , 
is significantly lower than in L1544 at radii 55 5,000 AU, 
which reflects the different density profile (see right bot- 
tom panel in Fig. [T5"|). 

The present data, together with previous work, show 
that there are significant chemical variations among ap- 
parently similar cloud cores and that CO is not always 
heavily depleted when the volume density becomes larger 
than a few xlO 4 cm -3 (as found in L1544, L1498 and 
L1517B; Caselli et al. 2002, Tafalla et al. 2002, 2004). 
Indeed, the case of L1521E, a Taurus starless core, where 
the centr al density is 10 5 cm" 3 but no CO freeze-out is 
observed (|Tafalla fc Santiagol2004bf ). suggests that cloud 
cores in similar dynamical stages can have different chem- 
ic al compositions. This point has been further discussed 
by I Lee et al.1 (|2003f ) , who underline the importance of the 
environment in setting the chemical/dynamical stage of 
a core, so that a core like L1521E (and L1689B) may 
have experienced a recent contraction phase, where the 
chemistry has not yet had the time to adjust to the new 
physical structure. On the other hand, the Bok Globule 
B68, being close to equilibrium, may have achieved the 
present structure a long time ago, so that both CO and 
N 2 H + had time to freeze-out (as found by iBergin et al.l 
pOOl L 

The present detailed study of TMC-1C adds a new 

9 The symbol Fd(CO), used here to indicate the CO deple- 
tion within the cloud, should not be confused with /d(CO), the 
observed (or integrated-along-the-line-of-site) CO depletion factor 
(/ D = J F u (l)dl; see also Crapsi et al. 2004). 
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piece to the puzzle: cores which are currently accreting 
material from the surrounding cloud appear chemically 
younger, with lower CO depletion factors. If the accret- 
ing cloud material, at densities ^10 4 cm -3 , is old enough 
(kJ 3 x 10 5 yr) to have formed observable abundances of 
N 2 H+ (as in TMC-1C), then N 2 H+ lines toward the dust 
peak will be bright. On the other hand, the chemical 
structure of cores such as L1521E (rich in CO, but poor 
in N-bearing species such as N 2 H + and NH 3 ) may be 
understood as young condensations which are accreting 
either lower density material (where the chemical times 
scales for N 2 H + formation are significantly longer) or 
material which spent only a small fraction of 10 5 yr at 
relatively high densities (~10 4 cm -3 ). The former hy- 
pothesis may be valid in environments less massive than 
those associated with TMC-1C (and it does not require 
large contraction speeds), whereas the latter hypothe- 
sis needs dynamical time scales shorter than or at most 
comparable to chemical time scales (which are about 10 4 
yr at densities of ~10 5 cm -3 , as can be found from the 
freeze-out time scale of species such as CO). 

In summary, TMC-1C, being more massive than L1544 
and other typical low-mass cores, has a larger reservoir of 
undepleted material at densities close to 10 4 cm~ 3 , where 
both CO and N 2 H + are abundant. Detailed chemical 
models suggest that TMC-1C must be at least 3 x 10 5 yr 
old, to reproduce the observed N 2 H + abundances across 
the cloud. Our simple chemical code tell us that the 
observed CO depletion factors can be reached in only 
10,000 yr. Therefore, the core nucleus is either signifi- 
cantly younger than the surrounding material or the sur- 
rounding (undepleted) material has accreted toward the 
core nucleus in the past 10,000 yr. In either case, this is 
evidence that the densest part of TMC-1C has recently 
accreted material. 

From the velocity gradients presented in Section 14.21 it 
is hard to see a clear pattern of flowing material towards 
the dust peak, but the "chaotic" pattern is reminiscent 
of a turbulent flow that may be funnelling towards the 
densest region, aided by gravity. In any case, the ex- 
tended inward motions deduced from N 2 H + observations 
(see Sect. I4.3[) is consistent with material at about (5- 
10)xl0 3 cm -3 currently accreting toward the dust peak 
position at velocities around 0.1 kms~ . It will be ex- 
tremely important to compare this velocity field with 
those predicted by turbulent simulations of molecular 
cloud evolution, especially considering the possibility of 
competitive accretion (Bonnell & Bate 2006). Observa- 
tion of CS will also be important to check our prediction 
of a larger "extended infall" velocity, when compared to 
L1544. 

7. SUMMARY 

A detailed observational study of the starless core 
TMC-1C, embedded in the Taurus molecular cloud, has 
been carried out with the IRAM 30m antenna. We have 
determined that TMC-1C is a relatively young core (t 
<; 3 x 10 5 yr), with evidence of material accreting toward 
the core nucleus (located at the dust emission peak). The 
core material at densities ^10 5 cm" 3 is embedded in a 
cloud condensation with total mass of about 14 M and 
average density of ~10 4 cm" 3 , where CO is mostly in 
the gas phase and N 2 H + had the time to reach the ob- 
served abundances of ~10 -10 . The overall structure is 



suggestive of ongoing inflow of material toward the cen- 
tral condensation. In addition, we have found that: 

1. N 2 H + (l-0) lines show signs of inward asymmetry over 
a region of about 7,000 AU in radius. This is the most ex- 
tended inward asymmetry observed in N 2 H + so far. The 
data are consistent with simple two-layer models, where 
the line-of-sight component of the relative (infall) veloc- 
ities range from ^0.15 kms" 1 (toward the dust peak) to 
~0.05 kms~ (at a distance from the dust peak of about 
7,000 AU). 

2. CO isotopologues and N 2 H + show increasing deple- 
tion as Ay increases and Td decreases. The amount of 
CO depletion that we observe is a factor of ~5 lower than 
that of L1544, whereas N 2 H + column densities are only 
a factor of two lower. Also, N 2 H + show clear signs of 
moderate depletion toward the dust peak position. 

3. The gas temperature determined from C 18 0(2-l) is 
12 K at the dust peak, indicating that CO is not trac- 
ing the dense (n(H 2 ) > 5 x 10 4 cm" 3 ) and cold (T kin < 
10 K) regions of dense cores. The C 18 0(2-l) excita- 
tion temperature drops outside the dust peak, and this 
is consistent with a roughly constant kinetic tempera- 
ture and a dropping volume density (traced by CO iso- 
topologues) from ~ 4 x 10 4 cm" 3 toward the dust peak 
to ~ 6 x 10 3 cm" 3 at a projected distance from the dust 
peak of about 11,000 AU (no high S/N data are available 
to probe larger size scales). 

4. N 2 H + (l-0) line widths are constant across the core, 
which is consistent with pre vious NH3 measurements 
(jBarranco fc Good man 1998), but different from what 
has been found with N 2 H + and N 2 D + observations 
of L1544 and L1521F (jCrapsi et al.l l2005h . where line 
widths are increasing toward the core center. The in- 
crease i n line width with radius seen in C 18 and N 2 H + 
in B68 (|Lada et al.ll2003D is not seen in TMC-1C, and 
unlike B68, the C 18 line width is significantly larger 
than the N 2 H+ line width throughout TMC-1C. This is 
consistent with the fact that TMC-1C, unlike B68, is em- 
bedded in a molecular cloud complex, so that CO lines 
trace more material along the line of sight of TMC-1C. 

5. The velocity field that we see in TMC-1C does not 
show the global signs of rotation that were seen in NH3 
observations over a somewhat differe nt area at arcminute 
resolution in I Goodman et aTl (|1993f) . Nevertheless, one 
portion of TMC-1C encompassing the dust peak posi- 
tion does have a more coherent velocity field, suggestive 
of solid body rotation with magnitude ~4 kms~ 1 pc" 1 , 
in the tracers C 17 O(l-0), C 17 0(2-l), C 18 0(2-l) and 
N 2 H+(l-0). 

6. The observed chemical structure of the TMC-1C core 
can be reproduced with a simple chemical model, adopt- 
ing the CO and N 2 binding energies recently measured 
in the laboratory. We argue here that "chemically young 
and physically evolved" cores like L1521E and L1698B 
(those with low CO depletion, faint N 2 H + lines, cen- 
tral densities above 10 5 cm -3 and centrally concentrated 
structure) have lower density envelopes than TMC-1C 
in which the N 2 H + abundance did not have the time to 
reach equilibrium values. On the other hand, "chem- 
ically and physically evolved cores" like L1544, L694- 
2 and L183 (those with high CO depletion and bright 
N 2 H + lines) are likely to have lower rates of accretion of 
material from the envelope to the nucleus than in TMC- 
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1C (or it has ended), and with a core nucleus undergoing 
contraction. Finally, "chemically evolved" but less cen- 
trally concentrated cores (e.g. L1498, L1512, B68), can 
just be older objects (age <; 10 6 yr), close to equilibrium, 
as suggested by Lada et al. (2003). In the case of TMC- 
1C, there is evidence that the core is at least 3 x 10 5 
yr old and has recently accreted less chemically evolved 
material. 

More comprehensive chemical models, taking into ac- 
count the accretion of chemically young material, as well 
as a comparison between the observed velocity patterns 
and turbulent models of cloud core formation are sorely 
needed to test our conclusions. 
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APPENDIX 
ERROR ESTIMATES ON T EX 



From the equation of radiative transfer (see eq. [9|), once T m b, r and the corresponding errors are known, the error 
on T cx (c>t ox ) can be determined following the rules of error propagation: 



2 / dT cx 

(Jrp = \ 7^; (Tx,,, 



nib 



3T ex 



where <7T mb an d o T are the errors associated with T m b and r, respectively. 
The expression of T cx is found by inverting eq. [9] 



T 



where 



A = 



Thus, the partial derivatives in eq. lAll are: 

dT cx 



]n[T /A + 1] 
a 2 b 



5T mb 

dT ex 
dr 
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where, 



a = T 
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c=J w (T h& ) 
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TABLE 1 

Observing Parameters of IRAM Spectra 



Transition 


T S y S 


Spectral Resolution 


FWHM 


B eff 


Frequency 


VLSR 




Kelvin 


km s~ 1 


arcseconds 




GHz 


km s — 1 


C 17 O(l-0) 


353 


0.052 


18.4 


0.66 


112.3592837 a 


5.2 


C 17 0(2-l) 


1893 


0.052 


9.2 


0.40 


224.7143850 a 


5.2 


C ls O(2-l) 


1050 


0.053 


9.4 


0.41 


219.5603541 a 


5.2 


C 34 S(2-1) 


202 


0.061 


21.4 


0.72 


96.4129495 a 


5.2 


DCO+(2-l) 


479 


0.041 


14.3 


0.54 


144.0773190 a 


5.2 


DCO+(3-2) 


622 


0.054 


9.5 


0.42 


216.1126045 a 


5.2 


N 2 H+(l-0) 


185 


0.031 


22.1 


0.73 


93.1737725 b 


5.2 



taken from lDore et al.l J2004D 



TABLE 2 

Fit Line Parameters at Dust Peak 



Transition 


VLSR 


/ Tdv a 


FWHM 


rms 




kms -1 


K kms" 1 


kms 


K 



C 17 O(l-0) 5.33 ± 0.03 0.37 ± 0.07 0.45 ± 0.11 0.089 

C 17 0(2-l) 5.49 ± 0.01 0.71 ± 0.05 0.43 ± 0.04 0.144 

C ls O(2-l) 5.20 ± 0.01 2.07 ± 0.01 0.41 ± 0.01 0.250 

C 34 S(2-1) 5.26 ± 0.04 0.17 ± 0.03 0.50 ± 0.10 0.092 

DCO+(2-l) 5.14 ± 0.01 0.28 ± 0.03 0.16 ± 0.02 0.204 

DCO+(3-2) 5.25 ± 0.02 0.34 ± 0.04 0.26 ± 0.04 0.179 

N 2 H+(l-0) 5.20 ± 0.01 0.20 ± 0.03 0.19 ± 0.03 0.104 

a Integrated intensity over entire spectrum (all hyperfine compo- 
nents included in the case of C 17 and N2H^). 



TABLE 3 

Fit results for the five spectra in Fig. [9] 



Offset 




TTOT 


Av 


Vf-Vt 


(arcsec) 


K 




kms -1 


km s 


-40,60 B a 


4.5 


17 


0.25 


0.0 


-40,60 F b 


3.3 


10 


0.15 


-0.05 


-20,40 B 


4.5 


15 


0.25 


0.0 


-20,40 F 


3.5 


15 


0.20 


0.05 


0,20 B 


4.4 


23 


0.17 


0.0 


0,20 F 


3.3 


10 


0.23 


0.15 


20,0 


4.4±0.3 


23±2 


0.17±0.01 




40,-20 


4.9±0.3 


8.8±0.8 


0.22±0.01 





a B=background layer F=foreground layer 



TABLE 4 

Molecular Transition Constants and Assumed Abundances 



Transition 


A a 


gb.c 


Si 


9u 








GHz 






n/n H2 


C 17 O(l-0) 


6.697E-8 


56.179990 


1 


3 


4.7 X 10~ 8c 


C 17 0(2-l) 


6.425E-7 


56.179990 


3 


5 


4.7 X 10" 8c 


C ls O(2-l) 


6.011E-7 


54.891420 


3 


5 


1.7 x 10~ 7f 


C 34 S(2-1) 


1.600E-5 


24.103548 


3 


5 


1.3 x 10- 10g 


DCO+(2-l) 


2.136E-4 


36.01976 


3 


5 


2.8 x lO" 1011 


DCO+(3-2) 


7.722E-4 


36.01976 


5 


7 


2.8 x lO" 1011 


N 2 H+(l-0) 


3.628E-5 


46.586867 


1 


3 


1.4 x 10~ lf,i 



a Einstein A coefficients are taken from the Leiden Atomic and Molecu- 
lar Database JSchoier et al.H2005f) b IFrerking et al.lll982T) c JGottlieb et al.l 

I2003H Standard molecular abundance, taken from literature except 
for N 2 H +e ICrapsi et all 120041) f l IGoldsmith et al.l [l99^ s l ITafalla et"aT1 

I200ai h I ILee et aUPOOa) ' Average value derived across TMC-1C. 
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TABLE 5 

Gas Column Density and Depletion at 
Dust Peak 



Transition 


N 


fo a 


Percent Error b 




cm -2 






C 17 O c 


1.0E15 


2.8 


5 


C 17 O d 


7.8E14 


3.6 


7 


c i8 d 


2.6E15 


3.8 


3 


C 34 s d 


5.4E11 


12 


20 


DCO+ d 


1.3E12 


10 


14 


DCO+ e 


2.8E12 


5.9 


13 


N 2 H+ C 


9.4E12 


0.9 


6 



a Depletion from dust- derived Nh 2 — 5.9 X 10 22 , 



N tot and X can Error from noise in spectrum, 
not including ~20% calibration uncertainty De- 
rived from (1-0) transition 1 ^ Derived from (2-1) 
transition Derived from (3-2) transition 
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Fig. 1. — The spectra taken at the (0,20) position, which is the peak of the dust emission map. The velocity scales on the spectra 
showing hypcrfine structure are only correct for the marked component, whereas the enlarged portion of the spectra has been recentered 
at the frequency of the isolated component. The vertical dotted line shows the LSR velocity of the weakest component of the N2H+(l-0) 
spectrum, as determined by a Gaussian fit, for comparison. 
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Fig. 2. — Spectra of the Fi,F = 2,3 — * 1,2 hyperfinc component of N2H+ across the whole TMC-1C map ped r egion. Red boxes mark 
spectra consistent with inward motions, whereas blue boxes mark spectra with outflow signatures (see Sect. l43t . Colored dots indicate 
those positions displayed in Fig. [9] The red-filled box mark the position of the continuum dust emission peak. The asymmetry in the high 
density tracer N2H+ is observed in an extended region, larger than that previously found in other evolved pre— stellar cores. See Fig. |4]for 
contours of column density based on dust emission. 
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Fig. 3. — Profile of the main Fi,F = 23 — > 12 hyperfine component of N2H+(1— 0) (white histogram) superposed on the spectrum of the 
weak FiF = 10 — * 11 component (yellow histogram) at the dust peak position. Note that the weak component is symmetric and peaks 
where the absorption of the main component is steepest, suggesting that this component is optically thin and can be used to estimate the 
column density and line width where the main group of hyperfines is affected by self-absorption. 





w 

K km/ s 

Fig. 4.— Integrated intensity maps of (top left) C 17 O(l-0), (top middle) C 17 0(2-l), (top right) C 18 0(2-l), (bottom left) DCO+(2-l) 
and (bottom right) N 2 H+(l-0). The C 17 O(l-0), C lr O(2-l) and DCO+(2-l) integrated intensities have been scaled up by factors of 1.5, 
2.0 and 5.0, respectively. For all spectra showing self-absorption, the N2H + integrated intensity has been calculated from the thinnest 
hypcrfine component. The contours show the visual extinction derived from the dust emission, at levels of Ay = 10, 25, 40, and 55 mag. 
The lo- uncertainties in the integrated intensity are ~0.05, ~0.06, ~0.10, ~0.07 and ~0.08 K km s^ 1 in C 17 O(l-0), C 17 0(2-l), C 18 0(2-l), 
DCO + (2— 1) and N2H + (1— 0) respectively. The positions with large N2H+ depletion and low column density have been marked with red 
crosses. 
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Fig. 5. — (Top Panel) N2H + (l-0) observed line width across the TMC-1C core. (Bottom Panel) N2H+(l-0) observed line width (in km 
s _1 ) vs. distance from the peak of the dust column density map. For those positions where the self— absorption is present, we used the 
width of the weak component as representative of the intrinsic line width of the (1-0) line. The crosses show the values at each position 
in the map with S/N > 3, and the large squares are averages of the data with 15" bins. The horizontal dashed line shows the thermal 
FWHM at 10 K for reference. 
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Fig. 6. — The non-thermal line width plotted against the thermal line width for each transition observed, at the position of the dust 
emission peak. The thermal line width is calculated from an assumed gas temperature of 10 K. The error in the non-thermal line width is 
derived using the error in the FWHM of the Gaussian fit to the line. The line withs of C 34 S(2-1) and DCO+(3-2) are highly uncertain, 
because of the stron g absorption observed (see Fig ,[2]l and the plotted errors do not take the absorption into account. The NH3 intrinsic line 
width is taken from Barranco & Goodman (1998). To make a fair comparison, all of the spectral line maps have been spatially smoothed 
to the ~1' resolution of the NH3 map. 
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Fig. 7.— Top Panels) The velocity field for C 17 O(l-0), C 17 0(2-l) and C 18 0(2-l). The black contour shows the A v = 10 mag boundary, 
and the arrows show the direction and magnitude of the line-of-sight velocity gradient, pointing from blue to red. The red box encloses 
a region with a velocity gradient similar to that of rotation as seen in the isotopologucs of CO. In each map the velocity shown is the 
difference between the measured VLSR at each position and the velocity at the (0,0) position. The lengths of the arrows are proportional 
to the magnitude of the gradient, and is scaled to the labeled arrow. {Bottom Panels) velocity field for DCO + (2— 1) and N2H + (1— 0) in 
TMC-1C. The red box encloses a region with a velocity gradient in iso topologues of CO that looks similar to that of rotation. The white 
box shows the region observed in NH3 by Barranco & Goodman (1998). 



22 



Schnee, S. et al. 



E 



1.5 l.i.5 0/5. 

km/s/pc 



D 



C 17 (2\ :) 




At 

C07(2 :) 



c 17 o ; 



C 18 (2-1) 



(1-0 



/ 

N H + (1-0)/ 



Fig. 8. — The overall velocity gradient for each tracer in the region bounded by —140 < RA j se t < 70 and —20 < DEC ff sc t < 70, the 
region in the red box in Fig. [7] and [7] The dotted circles show the magnitude of gradients equal to 0.5, 1.5, 2.5 and 3.5 km s _1 pc _1 . 
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Fig. 9. — N2H + (l-0) spectra (black histogram) along a South-East (bottom panel) to North— West (top panel) strip (see spectra labelled 
with colored dots in Fig- EJ . Filled red histograms are two-layer model results, whereas light blue filled histograms are simply fits to the 
line, assuming one emitting layer and constant T ex for the seven hyperfines. Offsets from Fig.[2]are in the top left of each panel. Numbers 
in the top right are the Vlsr velocities (in kms -1 ) from gaussian fits to the hfs structure. The parameters of the fits are in Table |31 Note 
that the redshifted absorption (and thus the infall velocity) is largest toward the dust peak and that the self-absorption is present towards 
the North- West but not towards the South— East. 
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Fig. 10. — Right ascension (left) and declination (right) cuts cross the TMC-1C core, passing through the dust peak position (black 
dotted line), of Ay (black curve), TV^^H" 1 ") (red curve), and X(N2H + ) (blue curve) normalized to their maximum values (see text). The 
blue dashed line is the abundance value observed across L1544 (see Vastel et al. 2006). X(N2H + ) shows an anticorrelation with the dust 
profile, suggestive of some N2H+ depletion toward the dust center. 
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Fig. 11. — The dust temperature plotted against the column density, derived from the dust emission maps at 450, 850 and 1200 /an. 
The spatial resolution of the data is 14", resolving the innermost region of the core, with Ay ~ 90 mag and T j UBt ~ 5 K. The solid red 
line shows the — Ay relation predicted for an externally heated pre-protostellar core by Zucconi ct al. (2001, Eq. 26). Departures from 
spherical symmetry are probably causing the observed discrepancy between observations and (spherically symmetric) model predictions. 
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Fig. 12. — {left) Excitation temperature of the C ls O(2— 1) line. If this line is in LTE, the present data suggest a gas temperature decrease 
from about 12 K at the dust peak, to about 10 K at > 1' away from the dust peak (equivalent to a projected distance of >8,000 AU). But 
we believe that this is caused by the density drop traced by the C ls O(2— 1) line (see text for details), (right) Brightness temperature ratio 
of the C 17 O(l-0) and C 17 0(2-l) lines as a function of distance from the dust (or Ay) peak. The ratio increases with distance, suggesting 
that some physical property (T^ n and/or n(H2)) changes as well. 
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Fig. 13.— Ratio of the C 17 O(l-0) and C 17 0(2-l) brightness temperatures as a function of ( Top) gas temperature, for a fixed value of 
the volume density (n(H2) = 2xl0 4 cm -3 ), and (Bottom) volume density, for a fixed value of the gas temperature (T^i n = 11 K). In both 
panels, 7V(C 17 0) = 10 15 cm -2 and At) = 0.4 kms" 1 . The dashed horizontal lines enclose the observed range of T m t> ratios (see Fig. 1121 
right panel). Note that observations are both consistent with a gas temperature decrease (and high constant density), as well as with a 
volume density decrease (and constant gas temperature) away from the dust peak. We believe that the latter is a more plausible solution 
than the former. 
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Fig. 14. — {left) The depletion factor derived from molecular transitions observed with the IRAM 30m plotted against the column density 
derived from dust emission maps at 450, 850 and 1200 /im. The median percent error in the depletion (calculated from the RMS of the 
spectrum) is given in each panel, (right) Same as left, except that the depletion factor is plotted against the dust temperature. Those 
points with anomalously high N2H+depletion factors (iVjy 2 < 3 X 10 22 and fr> > 0.6) are shown in red, and their positions are shown with 
red crosses in Fig. [4] The N2H+ depletion factor goes below 1 because of our arbitrary choice of the N2H+ column density (see text). 
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Fig. 15. — (Left panel) Fractional abundances of important molecular ions as a function of radius, calculated by the best— fit model 
for TMC— 1C (i.e. the one that best reproduces the observed column densities, see text). The red curves are depletion factors within 
TMC-1C (-Fb(CO), see the scale in the right y-axis) and, for comparison, within L1544 (dashed red curve, from Vastel et al. 2006). The 
discontinuity in the fractional abundance profiles i s due to the particular density profile used. (Right panel) Temperature and density 
profiles in TMC-1C, derived by Schncc ct al. (2007) and used in the chemical model. Dashed curves refer to L1544 (adopted by Vastel et 
al. 2006). 



